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Abstract 

The onset and bifurcation points of the n-cycles of a polynomial map are located through a 
characteristic equation connecting cyclic polynomials formed by periodic orbit points. The minimal 
polynomials of the critical parameters of the logistic, Henon, and cubic maps are obtained for n 
up to 13, 9, and 8, respectively. 



I. INTRODUCTION 



Consider the logistic map [TJ [2] : 

x k +i = f(x k ) =rx k (l- x k ). (1) 

If we iterate Eq. ([I]) from k — 0, what does the resulting sequence xq, x\ = f(xo), X2 = f{xi), 
. . . look like? We can visualize the sequence on the cobweb plot, see Fig. [I] for examples. 
Starting from (xo,xo) on the diagonal, each vertical arrow takes (x k , x/-) to (xk,y), where 
V — f( x k) = %k+i] the next horizontal arrow then reflects (xk,y) to (y, y) = (xi~+i,Xk+i), 
which starts the next iteration. 



( a ) Fixed point, r= 2.8 ( b ) 2-cycle, r= 3.4 ( c ) 3-cycle, r= 3.84 




FIG. 1. Cobweb plot of the logistic map. 



Three outcomes are possible: (i) a fixed point, which is a constant (including infinity), 
e.g., Fig. [l|a); (ii) a periodic cycle, which is a self-repeating pattern, e.g., Figs. [l](b)-(e); or 
(iii) a chaotic trajectory, e.g., Fig. []Jf). We will focus on the first two cases here. 

A fixed point is a solution of x* = f(x*). If x deviates slightly from x* and the sequence 
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still converges to x* , we call it stable. For a differentiable /, a stable fixed point requires 
\f'(x*)\ < 1 to reduce deviations in successive iterations pQ. 

In an n-cycle, n is the smallest positive integer that allows X\ = x n+ \ = f n (xi), where 
f n is the nth iterate of /, e.g., f 3 (x) = f(f(f(x))). Thus, any Xk in an n-cycle of / must 
be a fixed point of f n [the reverse is, however, untrue, for a fixed point of f n can also be 
a fixed point of f d as long as d\n: if f d (x) = x, then f n (x) = f d (- ■ ■ f d (x) ■ ■ ■ ) = x\ . We 
can therefore classify a cycle as stable or unstable by the corresponding fixed point of f n : a 
stable cycle requires | ^/"(^i) | < 1, or by the chain rule, 

f'(x l )...f(x n )\<l, (2) 

where the n points within the cycle, or the orbit. Further, the onset and 

bifurcation points are defined at the loci where ^/"(^i) reaches +1 and —1, respectively 

&■ 

The outcome of the iterated sequence of course depends on the parameter r. Below we 
will present an algorithm to identify all regions of r that allow stable n-cycles. 



II. LOGISTIC MAP 

We will illustrate the algorithm on the logistic map [D[2], defined in Eq. Q. If r is real, 
we will find windows (r , r&), within which stable n-cycles can exist. Here, if r > 0, then 
fa i r b) are the onset (bifurcation) points. There are generally multiple such windows even 
for a single n, but all onset points satisfy the same polynomial equation, and all bifurcation 
points another. Our goal is thus to find the two polynomials for a given n. 

To simplify the calculation, we first change variables [3] by 

x (ncw) ^ r ( x (old) _ R ^ r ( r _ 2 )/ 4; 

and rewrite the map, in terms of x (new) , as 

x k+1 = f(x k ) = R - x k 2 . (3) 

We will solve the cycle boundaries as zeros of the polynomials of R, and the corresponding 
polynomials for r can be obtained by R — > r(r — 2)/4. 
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A. Overall plan 



We solve the problem in two steps. Since the n-cycles form a subset of the fixed points 
of f n , we will first find the polynomials at the stability boundaries of the fixed points of f n 



(Sections II B to HE), then remove contributions from shorter <i-cycles (d\n) (Sections II F 



and [IIGJ). 

Let us consider the first step of finding the fixed points of /". At the first glance, the 
problem can be tackled by brute force: we can solve Eqs. ^ and express xi, . . . ,x n in terms 
of R, and then plug the solution into (j2|. The result contains R only (no Xk), and is therefore 
the answer. But since Eqs. ^ are nonlinear, it quickly becomes impossible for n > 2, as the 
degree of polynomials grows exponentially; thus it is nontrivial to reduce the final equation 
of R into a polynomial one. Nonetheless, on a computer, one can construct a Grobner basis 
[I] to automate the reduction. The approach, albeit straightforward, does not exploit the 
cyclic structure of Eqs. ([3]), can thus be improved by the following alternative. 

Instead of solving Eqs. (|3j) for Xk, we will derive a set of homogeneous linear equations of 
cyclic polynomials of x^ (an example of a cyclic polynomial is X1X2 + X2X3 + ■ ■ - + x n xi). Now 
the matrix formed by the coefficients of the homogeneous linear equations must have a zero 
determinant, for the cyclic polynomials are not zeros altogether. Thus, the zero-determinant 
condition gives the needed polynomial equations of R, whose roots contain all fixed points 
of f n . This completes the first step. 

For the second step, we show that short cycles serve as factors in the polynomials obtained 
above, and thus can be readily factored out. 



B. Cyclic polynomials 

A polynomial is cyclic if it is invariant under the cycling of variables x\ — > X2, x<2 — > 
X3, . . . , x n — > xi, e.g., a(x) = X1X2 + X2X3 + X3X4 + X4X1 and 6(x) = X1X3 + £2X4, for n = 4, 
where x = {xi, . . . , x n }. A cyclic polynomial should not to be confused with a symmetric 
polynomial, which is invariant under the exchange of any two Xk and Xj (j 7^ k); e.g., a(x) 
and 6(x) are not symmetric, but a(x) + 6(x) is. 

A cyclic polynomial can be generated by summing over distinct cyclic versions of a simpler 
polynomial of Xk, or a generator, e.g., Xix 2 is a generator of a(x); X\X 3 is that of 6(x); and 
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X1X2 + |xiX3 is that of a(x) + 6(x); note that the coefficient before £1X3 is 1 in the second 
case for there are only two distinct versions, but is | in the third case for there are four. 

We now consider cyclic polynomials generated from a monomial of unit coefficient, such 
as a(x) and 6(x), but not a(x) + 6(x). They can be systematically labeled as follows. We 
pick the monomial generator, which can be written as x\ ei X2 e2 ■ ■ ■ x n e ", then form a sequence 
p of indices with e\ l's, e 2 2's, . . . , e n n's; the corresponding cyclic polynomial is denoted by 

C p (x), e.g., Ci 2 (x) = xiX2 + X2X 3 H hx n Xi and C m (x) = Xi 2 x 2 + x 2 2 a;3H \-x n 2 Xi. We 

omit the length n in this notation, for we will mostly work with a fixed n at a time. Since a 
cyclic polynomial can have multiple generators, e.g., both x\x 2 and x 2 x 3 are generators of 
Ci 2 (x) (assuming n > 3), we pick the one that corresponds to the smallest p in the sense 
of lexicographic order, e.g., we choose x\x 2 instead of X2X3 for Ci 2 (x), because 12 < 23. 
Finally, we add Cq(x) = 1 for completeness. 

C. Square- free cyclic polynomials 

We further restrict ourselves to a subset of square-free cyclic polynomials, which have 
no square or higher powers of any e.g., Ci 2 (x) = x±x 2 + ■ ■ ■ is square- free, Cn 2 (x) = 
X\ 2 X2 + • • • is not, see Table [i] for more examples. Obviously, the label p of a square- 
free polynomial has no repeated index. We denote the set of all square-free p by B = 
{0, 1, 12, 13, ... , 123, 124, . . . , 12 . . . n} such that its size \B\ equals the number N& of square- 
free cyclic polynomials. 

We first show that the square-free cyclic polynomials serve as a basis for expanding cyclic 
polynomials: 

Theorem 1. For the logistic map Eq. |5p ; any cyclic polynomial K(x) formed by the n-cycle 
points x = {xi, . . . ,x n } is a linear combination of the square-free cyclic polynomials C p (x); 

where B = {0, 1, 12, 13, . . . , 12 . . . n} is the set of indices of all square-free cyclic polynomials, 
and f p (R) are polynomials of R. 

Proof. We show the theorem by the following square-free reduction. Given a cyclic polyno- 
mial K(x), we recursively apply Eq. Q as x\ — > R — Xk+i, until all squares or higher powers 
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TABLE I. Square- free cyclic polynomials C p (x) for the logistic map (n > 5). 



C P (x) 



Generator^ 



Necklace^ 



p as an index set* 




1 

12 
13 



1 

xi+ x 2 + ■ ■ 

X\X 2 + X 2 X 3 + 
X1X3 + X2X4 + 



Xr, 



x\ (or x 2 , ■■■) 
x n xi xix 2 (or x 2 x 3 , . . . ) 

x n x 2 xix 3 (or x 2 x±, . . . ) 



0...0 
10... 
110... 
1010... 



{1} (or {2},...) 
{1,2} (or {2,3}, ...) 
{1,3} (or {2,4}, ...) 



123 X1X2X3 + X2X3X4 + ■ ■ ■ + x n x\X2 X1X2X3 (or X2X3X4, . . . ) 1110... {1, 2, 3} (or {2, 3, 4}, 

124 X1X2X4, + X2X3X5 H + x„XiX3 X1X2X4 (or X2X3X5, . . . ) 11010... {1, 2, 4} (or {2, 3, 4}, 



12... n 



X\X 2 ■ ■ ■ x r 



X1X2 ■ ■ ■ x n 



11 ... 1 



{l,2,...,n} 



t Alternative generators are shown in parentheses. 

* The corresponding binary necklaces. 

* The label p of a square- free cyclic polynomial has no repeated indices; so the indices can be cast to a set. 



of Xk are eliminated. The process will not last indefinitely for each substitution reduces the 
degree in Xk (no matter which k) by one. Since the original polynomial is cyclic, so is the 
reduced one. All terms that involve no x^ are collected to serve as the coefficient before 
Co(x), which is 1. Since no square or higher powers of Xk can survive the reduction, all 
cyclic polynomials in the final result are square-free. The coefficients are polynomials of R, 
for R is the only variable introduced by the substitutions. □ 

For example, for n = 2, the cyclic polynomial K(x) = x\ 2 x<i + X2 2 X\ can be written as 
K{x) = (R + l) Ci(x) — 2RC (x) for x 2 x 2 = Rx 2 — x 2 2 = Rx 2 — R + 21 and x 2 x\ = 
Rx\ — R + x 2 . 

Theorem [T] shows that any cyclic polynomial can be expanded as a combination of the 
square-free ones, which serve as a basis. Below we show that at the onset and bifurcation 
points, the square-free cyclic polynomials C p (x) are themselves linearly connected by an 
N$ x Nb matrix equation. The determinant of matrix must vanish, and this condition yields 
the solution of the problem. 
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D. Algorithm for locating fixed points of / 



We first observe that the derivative of f n is a cyclic polynomial: 

A(x) = ^/ n (zi) = f'M . . . fix,) = (-2) n Xl ...x n . (4) 

Now, for any p, A(x) C p (x) is also a cyclic polynomial, since the product of two cyclic 
polynomials is cyclic too. We can therefore expand it by Theorem [T] as 

A(x)C p (x) = ^T P9 ( J R)C (? (x), (5) 

where T pq (R) is a polynomial of R, and p,q G B. 

By Eq. (|2]), at the onset or bifurcation point, A(x) is equal to a number A = +1 or — 1, 
respectively; so Eq. ^ becomes a homogeneous linear equation of C p (x): 

AC p (x) = ^T M (E)C 9 (x), (6) 

or in matrix form, 

[AI-T(fl)]C = 0, @) 

where I is the N& x N& identity matrix, T(R) = {T pq (R)} is an Ng x Ng matrix, and 
C = {C p (x)} is an A^-dimensional column vector. 

Since a set of homogeneous linear equations has a non-trivial solution only if the deter- 
minant of the coefficient matrix is zero, we have 

A n (R,X) = |AI-T(#)| =0. (7) 

Here we have defined A n (R, A) as a polynomial of R and A, and we have also attached the 
subscript n, for later use with Ad(R,X), where d are divisors of n. Eq. ^ is a necessary 
condition since C p (x) cannot vanish altogether; and since it involves R only, the polynomial 
expansion of the determinant gives the answer to our problem. 
To summarize, we have 

Theorem 2. At the onset and bifurcation points, the square-free cyclic polynomials C p (x) 
are linear related by Eq. M) ; with A being +1 and —1, respectively, andT pq (R) the coefficients 
from the square-free reduction o/A(x) C p (x) with A(x) specified by Eq. Q). Thus, R at the 
two points are the roots of the polynomials A n (R,X = ±1) obtained from the characteristic 
equation Eq. (fTl). 



Remark 1. For complex R and x, A should be generalized from ±1 to any A = exp(z0), 
where (f) G [0, 2ir], the algorithm still applies. By increasing <p from to 2ir, we can trace 
a two-dimensional region of a complex R for stable cycles. These regions are bulbs in the 
Mandelbrot set, see ref. [5] and Fig. [2j 

Further, with A = 0, the algorithm determines the superstable point, at which the de- 
viation from a cycle point vanishes to the linear order after n iterations of /. But we 
have a better algorithm in this case: since at least one of the Xk is zero by Eq. Q, then 
f n (%k = 0) = x n+ k = provides the needed polynomial equation of R [I]. 



E. Examples 



We illustrate the above algorithm by cases of small n. It is still helpful to have a mathe- 
matical software verify some steps (e.g., in computing the determinants and their factoriza- 
tion). 

For n = 1, we have two square-free cyclic polynomials Co(x) = 1 and Ci(x) = x%) and 
A(x) = — 2x\ [we shall drop "(x)" below for convenience]. Thus A Co = — 2x\ = —2C\, 
ACi = -2x\ = -2R + 2 Xl = -2RC + 2C U or 









-4R-2A + A 2 , 



1 ) 



and Eq. reads 

A 2 
2R A - 2 

which is just the equation for a fixed point. The fixed point begins atA = +lori? a = — 1/4, 
and becomes impossible when A = — 1, or R b = 3/4. 

For n = 2, the cyclic variables are Cq = 1, G\ = x\ + X2, C12 = %i%2, an d A = 4x±X2- 
Thus, AC = 4C 12 , ACi = 4xi 2 x 2 + 4x x x 2 2 = 4R(xi + x 2 ) -4(a;i 2 + X2 2 ) = AR(x l +x 2 )- 
8R+A( Xl +x 2 ) =A(R+1)C 1 -8RC ,AC 12 = 4(R-x 2 ) (R-xx) = AR 2 C -ARd+AC 12 , 
or 



A 



( C \ 



4 

-BR A(R+ 1) 
4i? 2 -4R 4 



Co 
Ci 

v c 12/ 



8 



and Eq. ^ reads 



A -4 

8R A-4(i?+l) 
-4R 2 AR A - 4 

(4#-4 + \)[{AR - A) 2 -4A]. 



We only use the first factor (the choice will be explained later, same for the following cases). 
Setting it to zero yields R = 1 — A/4; A = +1 gives the onset value R a = 3/4 (r = 3) while 
A = — 1 gives the bifurcation value R& = 5/4 (r^ = 1 + v6)- Note that the onset point of 
the only 2-cycle is located at R — 3/4, where the fixed point bifurcates p] [compare Figs, 
[if a) and (b)] . 

For n = 3 (21 EHS], we have C = 1, C\ = X\ + x 2 + x 3 , G\ 2 = X\X 2 + x 2 x 3 + x 3 Xi, 
C123 = xix 2 x 3 , and A = —8x1X2X3. The square-free reduction yields 



A 



/ Co \ 

Ci 

\ c 123 J 







/ 

-2AR 8(72 + 1) -8/2 

24i? 2 -8R(R + 2) 8(R+1) 

\ -8i? 3 8i? 2 -8R 8 



-8\ ( C \ 

Ci 



/ 



Cl2 
V ^123 / 



(8) 



and Eq. Q reads: 



AO 8 

24# A-8(i?+l) 8i? 

-24i? 2 8R(R + 2) \-8{R + l) 

8i? 3 -8i? 2 8i? A- 



[64R 3 - 128R 2 - 8(A - 8)R - (A - 8) 2 ] (A 2 - 8A - 24RX - 64R 3 ) . 



Using the first factor, we find at the onset point A = 1, (R — |) (i? 2 — \R+ ^) =0 and 
its only real solution is i? a = 7/4 (r a = 1 + v8). At the bifurcation point A = — 1, the 
equation R 3 - 2R 2 + |i? - § = yields # 6 = ± ( | + ^ff - fv^Ol + + fv 7 ^ 



whose corresponding r = 1 + a/1 + 4i? is identical to that in ref. [3, IE]- 

For n = 4 [5], the cyclic variables are Co = 1, G\ = x\ + x 2 + x 3 + X4, C12 = X1X2 + x 2 x 3 + 
x 3 x 4 + x A x 1: C 13 = Xix 3 + x 2 x 4 , C123 = X1X2X3 + X2X3X4 + X3X4X1+X4X1X2 C1234 = £1X2X3X4, 
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and A = I6X1X2X3X4. Eq. ^ reads 



AO -16 

647? A -16(72 + 1) 167? -167? 

-647? 2 167?(7? + 2) A-16(7? 2 + l) -327? 16R 

-327? 2 167?(7?+1) -167? A - 16(7? 2 + 1) 

647? 3 -167? 2 (7? + 3) 167?(7? + 2) 327?(7? + 1) A - 16(72 + 1) 

— 16i? 4 167? 3 -167? 2 -167? 2 1672 A - 16 

= [40967? 6 - 122887? 5 + 256(A + 48) (7? 4 - 7? 3 ) - 16(A + 32) (A - 16)7? 2 - (A - 16) 3 ] 

[16(72 - l) 2 - A] [(16i? 2 + A) 2 - 16(27? + 1) 2 A] , 04) 

From the first factor, we have (47? - 5) [(47? + l) 2 + 4] [(47? - 3) 3 - 108] = at the onset 
point A = 1. It has two real roots: 7?^ = 5/4 (r' a = 1 + v^6 ~ 3.4495) for the cycle 
from period-doubling the 2-cycle [compare Figs. lVb) and (d)] , and 7? a = (3 + v / 108)/4 
(r a = 1 + ^4+ ^108 w 3.9601) for an original cycle [Fig. [lj(e)] . At the bifurcation 
point, A = -1, and 40967? 6 - 122887? 5 + 12032(7? 4 - 7? 3 ) + 8432T? 2 + 4913 = 0, which 
upon 7? — > r(r — 2)/4 yields the same polynomial obtained previously jH [TUHT2"] . The only 
two positive roots « 3.9608 and r' b w 3.5441 correspond to r a and r' a respectively. As a 
verification, the polynomials are alternatively derived in Appendix [A] 

The algorithm was coded into a Mathematica program, which was used to compute the 
polynomials for n up to 13. The polynomials for a general A and those at A = ±1 (onset 
and bifurcation points) are listed in Table |n] and Table IV respectively, for some small n. 



For complex 7?, A, and x, the method can also compute the region of stability for 7?, with A 
being exp(i0) [0 G (0, 27r)] instead of ±1; the results are shown in Fig. [2] for n up to 8. For 
polynomials of larger n, see the website in Section |Vj The representative r values are listed 
in Table [TTT1 



F. Minimal polynomial for the n-cycles 



The factors ignored in Section II E come from shorter <i-cycles whose periods d divide n, 
because A n (R, A), from the characteristic equation Eq. Q, is derived for all fixed points of 
and thus encompasses the shorter cycles as well. We filter the contributions from the 
shorter cycles by the following theorem. 
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TABLE II. Characteristic polynomials A n (R, A) of the fixed points of f n of the simplified logistic 
map Eq.©. 



A n (R, \ = 2 n X)/2 



nN B (n) f 



X 2 -X-R 




[iB^xj 2 -^] (R + X-l) 

: ^R 5 '] [(X - l) 2 +XR- R(R - l) 2 ] 
4 [iR^+f-^X] [(^JrJC^ipt+lfX] [R 6 - 3R 5 + {X + 3)(i? 4 - R 3 ) -(X + 2)(X- 1)R 2 - {X - l) 3 ] 

[X^JiXjM^r-R^ [ - R 15 + 8R 14 - 2SR 13 + (X + 60)i? 12 - (7X + 94)i? n + {3X 2 + 20X 

+116)i? 10 - (11X 2 + 33X + IU)R 9 + 2(3X 2 + 20X + 47)i? 8 - (2X 3 - 20X 2 + 37X + 69)i? 7 

5 

+ (3X - 11)(3X 2 - 3X + U)R 6 - (X - l)(3X 3 + 20X 2 - 33X - 26)R 5 + (X - 1) 2 (3X 2 + 27X + 14)i? 4 
+ (X - l) 3 (6X + 5)R 3 + (X - 1) 4 (X + 2)R 2 + (X - l) 5 R + {X - if] 

[iR^f^F^] \(R 3 -jQ^-^iS^nyx] \{R(R - ^y^X^X) 2 ^^ - 2) 2 X\ [R 27 - 13R 26 + 78R 25 
6 

+ (X - 293)i? 24 + ■ • • + (X - lf(X 2 + 10X + 3)R 3 + (X - 1) 7 (X + l)i? 2 - (X - l) 8 R + (X - if] 
[-R 7 -JXM^r^tf^^ [ - R 63 + 32R 62 - 496i? 61 + 4976i? 60 + {X - 36440)i? 59 - (30X 
208336)i? 58 + • • • + (X - 1) 15 (X 2 + UX + 5)R 3 + 2(X - l) 16 (X + l)R 2 + (X - 1) 17 R + (X - l) 18 ] 



Nein) = {l/n) J2d\n <t>{ n /d)2 d [Eq. (11 |] is the number of the square-free cyclic polynomials. The change 



of variable A — > X and the division by 2 nNB ^ make the polynomials more compact. The irrelevant factors 



from shorter cycles (see Section II F) arc struck out. The polynomials of r for the original logistic map 
Eq. jlj) can be obtained by R — » r(r — 2)/4. 



Theorem 3. The minimal polynomial P n (R, A) of all n-cycles is a factor of A n (R, A) [defined 
in Eq. ffi], and can be computed as 



P n {R,X) = H B d)C (R,X) 



(9) 



cd=n 



where Bd, c (R,X) = Ylk=i Ad(R, e 2km l c \ l l c ) , \ l ' c is a complex cth root of X, and /i(c) is the 
Mobius function. 
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TABLE III. Smallest positive r at the onset and bifurcation points of the n-cycles of the logistic 



map. 





Onset* 


Bifurcation* 


#* 


rat 


Onset* 


Bifurcation* 


#* 


1 


l.OOOOOOOOOOi 


3.0000000000i 


1 


8' 


3.9607686524 6 


3.9610986335i2 O 


1 


2' 


3.0000000000i 


3.4494897428i 


1 


8"' 


3.5440903596 6 


3.5644072661i2o 


1 


3 


3.8284271247i 


3.8414990075 3 


1 


9 


3.6871968733 240 


3.6872742105 252 


28 


4 


3.9601018827 3 


3.9607686524 6 


1 


10 


3.6052080669 472 


3.6059169323 495 


48 


4" 


3.4494897428i 


3.5440903596 6 


1 


10' 


3.7411207566i5 


3.7425706462 495 


3 


5 


3.7381723753n 


3.7411207566i5 


3 


11 


3.6817160194 10 i 3 


3.6817266457io2 3 


93 


G 


3.6265531617 20 


3.6303887000 27 


4 


12 


3.5820230011i959 


3.58 28117795 20 io 


165 


C 


3.8414990075 3 


3.8476106612 27 


1 


12' 


3.6303887000 27 


3.6321857392 20 io 


4 


7 


3.7016407642 57 


3.7021549282 63 


9 


12" 


3.8476106612 27 


3.8490363152 20 io 


1 


8 


3.6621089132io8 


3.6624407072i2 


14 


13 


3.6797024578 40 83 


3.6797038498 40 95 


315 



, or "' means a cycle undergoing the first, second, or third successive period-doubling, respectively. 



* The subscripts are the degrees of the corresponding minimal polynomial of R — r(r — 2)/4. 

* The number of similar cycles. 



Remark 1. The Mobius function fi(n) is (— l) fe if n is the product of k distinct primes, or 
if n is divisible by a square of a prime, //(n) = 1, —1, —1, 0, —1, 1, . . . , starting from n = 1. 
The fi(n) is useful for inversion: g{n) = 2~^d|n M 71 /*^) M^) ^ anc ^ om y ^ h{n) = Yld\n9^) 
PJ. 

Remark 2. B^ C (R, A) is a polynomial of A. Despite the argument A 1//c , the product 
U c k=1 A d (R,e 2km/c \ 1/c ) is free from radicals of powers of A 1//c , for it is invariant under A — > 
e 2 ™A; and deg A B d>c (R, A) = deg A A d (R, A). Particularly, B n>1 (R, A) = A n (R, A). 



Let us see some examples. For n = 1, there is no irrelevant factor in Eq. (7-1) and 
Px(R, A) = Bx,i(R, A) = Ai(R, A) = A 2 - 2A - 4R. 



For n = 2, since B 1>2 {R, A) = (A + 2y/\ - AR) (A - 2y/\ - 4R) = (4.R - A) 2 - 4A, 
P 2 (R, A) = A 2 (R, A) B 1>2 {R, A)" 1 = AR - 4 + A. 

For n = 3, one can verify that B h3 (R, A) = flLi A ^ e 2kni / 3 = A 2 - 24i?A - 8A - 
64.R 3 . So P 3 (R, A) = A 3 (R, A) B lt3 (R, X)' 1 = - [GAR 3 - 128R 2 - 8(A - 8)R - (A - 8) 2 ] . 
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Forn = 4, we have 16(i?-l) 2 -A = (4R-4+VA)(4R-4- and (16i? 2 + A) 2 -16(2i?+ 
1) 2 A = [(4i2 - v^) 2 - 4^] [(4R + yf\) 2 + 4VX] ■ Thus, the last two factors of Eq. (frQ 



can be written as B 2>2 (R, A) = nLi M{R, e kni y/\), and P 4 (R, A) = A 4 (i2, A) B 2>2 (R, A)" 1 = 
4096 J R 6 -12288i? 5 + 256(A + 48)( J R 4 -i? 3 )-16(A + 32)(A-16) J R 2 -(A-16) 3 . Note, B 1A (R,X) 
is unused for /i(4) = 0. 

The irrelevant factors for n up to 7 are listed in Table [TTJ 




FIG. 2. Stable regions of the n-cycles of the simplified logistic map Eq. Q with a complex R; 
obtained from P n (R, e**) = with <p G [0, 2vr] cf. Fig. 3 in [5]. 



Theorem[3]is not always necessary. For n > 4, P n (R, A) is readily recognized as the factor 



of A n (i£, A) with the highest degree in R, see Table [H| and Section III It can be, however, 
problematic, if P n {R, A) is solved for A = 1 instead of a general A, for P n (R, 1) itself can be 



further factorized, see Table IV Due to the technical nature of the derivation and subsequent 
discussions, the reader may wish to skip the rest of Section [IT] on first reading. 
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TABLE IV. Onset and bifurcation polynomials of the n-cycles of the simplified logistic map Eq. O 



Onset P n (R, +1) t.t 



Bifurcation P n (R, — 1) * 



-R 4 - 1 

(£^5) [tB^H^F^ [{Ri - 3) 3 - 108] 
-416i? 4 9 - 3404i? 4 8 + 20548i? 4 7 - 98258i? 4 6 



+370146i? 4 5 - 1171676i?4 4 + 3301996i? 4 3 



-7507332i? 4 2 + 15699857i? 4 - 28629151) 

r+3)(i?4 3 
3063651608241) 

5T 2 + 97i? 4 + 127) 
58165204 . . . 8504447) 



[Ri — Rn° + i?. 4 




-i? 4 + 3 
R 4 - 5 



-i? 4 3 + 8i?4 2 — 18i?4 + 81 



(i? 4 57 - 127i? 4 56 + 



-R 4 6 - 12i? 4 5 + 47i? 4 4 - 188i? 4 3 + 527i? 4 2 + 4913 



-.R4 15 + 32i? 4 14 - 448i? 4 13 + 3838i? 4 12 - 24008i? 4 n 



-118147i? 4 10 - 462764i? 4 9 + 1519712i? 4 8 



-4444924i? 4 7 + 11351480^4 6 - 26978787i? 4 5 
+58697100i? 4 4 - 88548768i? 4 3 + 149426046i? 4 2 
-313083144^ + 1291467969 



R 4 27 - 52i? 4 26 + 1248i? 4 25 - 18753i? 4 24 



5098317006250000i? 4 - 20711912837890625 



i? 4 63 + 128i? 4 62 - 7936i? 4 61 + 318464i? 4 60 



2427583 . . . 6441888i? 4 + 9786215 . . . 0031361 




/n 12 nil 
- 1 1L4 — -fl4 



18433i? 4 + 8191) 



13 (^4083 _ 8191jR4 4082 + 33529856^*08! 



-30826655683291995 . . . 27828275886475354111 



-i? 4 4095 + 8192i? 4 4094 - 33538048i? 4 4093 



7361199006999 . . . 96207964555264i? 4 

-29448390363448812 . . . 352154556569141249 



t i?4 = 4R. The polynomials of r for the original logistic map Eq. ( 1 1 can be obtained by R4 —> r(r — 2) 
* Factors for the n-cycles born out of shorter cycles (see Section II Hi are struck out. 



G. Counting cycles 

To show Theorem |3j we first find the degrees in A of A n (R, A) (Theorem [1]) and P n (R, A) 
(Theorem |5]) . By comparing the degrees, we then show that each Pd(R, A) with (d\n), after 
some transformation, contributes one polynomial factor to A n (R, A) (Theorem [6]) , and the 
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inversion of the relation yields Theorem [3j 



1. Number of the square- free cyclic polynomials 

To count the square-free cyclic polynomials, we establish a one-to-one mapping between 
the square-free cyclic polynomials and the binary necklaces (defined below). The task is 
then to count the latter. 

A binary necklace is a nonequivalent binary 0-1 string. Two strings are equivalent if they 
differ only by a circular shift. For example, for n = 3, there are 2 3 = 8 binary strings, 
but only four necklaces: 000, 001, 011 and 111, since 010 and 100 are equivalent to 001, so 
are 110 and 101 to 011. The period of a necklace, or a binary string, is the length of the 
shortest non-repeating sub-sequence, e.g., the periods of 1111, 0101 and 0001 are 1, 2, and 
4, respectively. Obviously, the period m divides n; and a period-m necklace encompasses m 
binary strings differed by circular shifts, e.g., 0101 represents both 0101 and 1010. 



TABLE V. Two ways of counting the iV(4) = 6 necklaces for n = 4. 





m = 1 


m = 2 


m = A 


2 d 


4>(n/d) 


<j){n/d)2 d 


d = 1 


X 4 1x4 






2 = 2 


2 


4 


d = 2 


00x2 Hx2 


10 x2 (01x2) 




2 + 2 = 4 


1 


4 


d = 4 


0000 1111 


1010(0101) 


1000(0100,0010,0001) 
1100(1001,0110,0011) 
1110(1101, 1011,0111) 


2 + 2 + 12 
= 16 


1 


16 


m Em|,ld|n ( M§) 


1- (2 + 1 + 1) =4 


2- (1 + 1) = 4 


4-1 = 4 




1 


N(n) ■ n 


2-4 


1-4 


3-4 


— > 


6 • 4 = 24 



Bold strings are necklaces; others are their cyclic versions. 

The subscript of a binary string means the number of repeats; e.g., X 4 means repeated four times, or 
0000; 10 X 2 means 10 repeated twice, or 1010, etc. 



To compute the number of the binary necklaces N(n), we construct a sum for the length- 
n binary strings and count it in two ways. Table |Vj shows the example for the n = 4 case. 
For each divisor d of n, we collect all length-n binary strings whose periods m divide d. 
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The total is 2 , for we have enumerated all binary strings whose periods divide d. We then 
weight them by the Euler's totient function <p{n/d). Here, <f)(m) gives the number of integers 
from 1 to m that are coprime to m, e.g., 0(1) = 1, 0(3) = 2 for 1 and 2, 0(6) = 2 for 1 and 
5. We repeat the process over other divisors d of n, and the resulting sum is Yld\n ( f ) ( n /^')^ d - 



We can count the above sum in another way. We recall that a period-m necklace always 
contributes m strings in the above process for a fixed d, and it does so for all multiples d of 
m. So the total weighted contribution by this necklace is 



where we have used the identity ^2 d /\ m i 4>{d') = w! , with d! = n/d and m' = n/m. Summing 
over the necklaces yields N{n)n. Thus, N(n) = (l-/n)J2d\n ( f ) ( n /d)'2 d - The process for a 



Back to our problem, there is a correspondence between the binary necklaces and the 
square-free cyclic polynomials. For a square-free cyclic polynomial, we construct a binary 
string according to its generator: if it contains x^, the kth character from the left is 1, 
otherwise 0. The resulting string corresponds to a unique necklace; the alternative generators 
give the circularly-shifted binary strings. The mapping is reversible, or one-to-one. For 
example, if n — 3, the square-free polynomials for 000, 100, 110 and 111 are Co = 1 
(generator: 1), C\ = Xi + x 2 + x 3 (generator: x±), C\ 2 = X\X 2 + x 2 x 3 + x 3 Xi (generator: 
Xix 2 ) and C123 = X\X 2 x% (generator: x\X 2 x^), respectively. Table [T] shows a few more 
examples. Thus, 

Theorem 4. The number N^n) of the square-free cyclic polynomials C p (x) (p G B) formed 
byx= {xi, . . . ,x n } is 



The process for a fixed d is exemplified by a row in Table |Vj 




(10) 



fixed necklace is exemplified by a column in Table |V} 




(11) 



d\n 



which is also equal to deg A A n (R, A). 



N(n) = 2, 3, 4, 6, 8, 14, 20, 36, 60, 108, 188, 352, 632, . . . , starting from n = 1. 
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2. Number of the n-cycles 



Theorem 5. The degree in A of the minimal polynomial P n (R, A) of all n-cycles is equal to 
the number of the n-cycles, and is given by 



Proof. Except a few special values of R, the iterated map f n generally has 2 n distinct 
complex fixed points, for otherwise f n (x) — x would have a repeated zero at any R, but at 
R = 0, —x 2 " — x has no repeated root; a contradiction. 

Each fixed point can be assigned to a point in a d-cycle, with d being a divisor of n. 
The assignment is both complete (for a <i-cycle point must also be a fixed point of f n ) 
and non-redundant (for there is no repeated fixed point of f n ). Since each of the L(d) 
d-cycles contributes d fixed points, we have 2 n = Yld\n L(d) d. The Mobius inversion yields 
L(n) = (1/n) Yld\n fi{n/d)2 d . This formula was known to several authors [141 IT5] . 

We now define A c as the value of A(x), evaluated at the cycle points x( c ) = {x^ c \ . . . , x n ^} 
of cycle c. Of course, A c is a function of R. If we assume that A c are distinct (see Remark 
1 below), then the minimal polynomial P n (R, A), as a polynomial of A, takes the form of 
ric(^ — ^c)- Thus, the degree of P n (R, A) in A must be the same as the number of the 
n-cycles. □ 

Remark 1. Although all trajectory points x^ c ^ are distinct in different cycles, the value 
A c of the polynomial A(x^ c ^) may happen to be the same. In this case, we shall find another 
cyclic polynomial Y(x) that has different values in the two cycles [Y(x) exists for otherwise 
x are the same in the two cycles], use A'(x) = A(x) + eF(x) to list Eqs. ([5]), then take the 
limit e — > 0. 

Remark 2. L(n) is also the number of aperiodic (i.e., period equal to n) necklaces of 
length n, e.g., for n = 4, out of the six necklaces, 0001, 0011, and 0111 are aperiodic; but 
0000, 1111, and 0101 are periodic. Since a period-c? (d\n) necklace is just an aperiodic one 
of length-d repeated n/d times, and each contributes d binary strings, we can count the 2™ 
binary strings of length n as 2 n = Yld\n L{d) d. The Mobius inversion leads to the same 




(12) 



d\n 



result of Eq. (12). 



L(n) = 2, 1, 2, 3, 6, 9, 18, 30, 56, 99, 186, 335, 630, . . . , starting from n = 1. 
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Since the period d of a length-n necklace always divides n, we have 



N(n) = J2m- (13) 

d\n 



We can also show this by explicit computation: 



^L(rf) = ^(l/ C /)^M^/c)2 c = (l/n)^2 c ^(^V 

d\n d\n c\d c\n d'\(n/c) 

= (l/n)^2 c 0(n/c) = iV(n), 

c I n 

where we have used 0(m) = X^d'im "which i s the inversion of m = $^ d /i m </>(rf')j an d 

Eq. (11 ). We will use Eq. (fl~3|) in proving the next theorem. 



3. Relation between P n (R, A) and A n (R, A) 

Theorem 6. The minimal polynomials Pd{R, A) of all d-cycles of periods d\n and A n (R, A) 
[defined in Eq. Pity] are related by 

A n (R,X)= ]jQ d7C (R,X), (14) 

cd=n 

where Q d ,c(R, X) = Ylk=i Rd(R, e 2fc7rt / c A 1//c ) is a polynomial of degree L{d) in A representing 
contributions from d-cycles. 

Proof. Since A n (R,X) represent all <i-cycles with d\n, and each cycle holds a distinct A, 
deg A A n (R, A) is at least ^2 d \ n L(d) according to Eq. (12), which is equal to N(n) by Eq. (13). 



Since deg A A n (R, A) = N(n) by Eq. (11), each n-cycle occurs exactly once in A n (R, A). 

In a d-cycle, we have A({xi, . . . , x n }) = (— 2) n xi . . . x n = A/, where A d ({xi, . . . , Xd\) = 
(— 2) d x\ . . . Xd and c = n/d. So the (i-cycle satisfies a polynomial Pd(R, A 1//c ) = 0, where 
A = A({xi, . . . ,x n }). This is, however, not a polynomial equation, and the radical A 1//c can 
be removed by the product Q d>c {R, A) = Ul=i p d(R, e 2kni/c X 1/c ) = 0. Now Q d>c {R,X) is a 
polynomial of A for it is invariant under A — > e 2m X, and thus free from radicals of the form 
A'/ c [if {I, c) ^ c] . And since deg A Q d ,c{R, X) = deg A P d (R, A) = it is also a polynomial 
of the lowest possible degree in A. 

Therefore, the product n c d=n Qd,c(R, A) can differ from A n (R, A) only by a multiple. Since 
Qn t i(R, A) = P n (R, A), and the coefficient of highest power of A is always unity in A n (R, A) 
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[see the definition Eq. (jTb] , we know by induction that the coefficients of the highest power 



of A in all P n (R, A) and Qd,c(R, A) are also unities. So the multiple is one, hence Eq. (14). □ 
We can now prove Theorem [3] as a corollary of Theorem [6j 

m 

B n , m (R,\) = n^n^e^^A 1 /" 1 ) 



m c 



= n n n p »d R > e^+^/^A 1 ^) 

c\n 1=1 k=l 

mc 

= n n p n/ c{R,^ /(me) x i/ime) ) = i[Q n /cmc{R, a). 

c\n k'=l c\n 

Taking the logarithm (formally) yields log-B njm = Yld\n 1°S Qd,mn/d, where d = n/c. The 
inversion is logQ njm = Y,d\n^( n / d ) ^°S B d,mn/d, or 

Qn,m(R, A) = Y\ B d)mn / d (R, \)^ n / d \ 
d\n 

which is reduced to Eq. ^ with m = 1 for Q ni i(R, A) = P n (R, A). 

H. Intersection of cycles and further factorization at the onset point 



Table IV shows that the onset polynomial P n (R, A = 1) for the n-cycles can be further 
factorized. This is because the intersection of an n-cycle and a shorter d-cycle (d\n, d < n) 
forces the two to share orbits (this cannot happen if d / n, for the orbits would be out of 
phase). Consequently, upon the intersection, P n (R, A) from the n-cycle has to accommodate 
Pd{R, A') from the ci-cycle, with A' being a primitive (n/d) th root of A. 

At the intersection, the shorter d-cycle is branched or "bifurcated" by (ra/d)-fold to the 
n-cycle. The simplest example is the first bifurcation point at R = 3/4 for d — l,n = 2, 



where the fixed point Eq. (7-1) bifurcates to the 2-cycle Eq. (7-2). The second bifurcation 



point at R = 5/4 for d = 2, n = 4 is similar, cf. Section II E 

We will show below that such branching generally can only happen at the onset of the 
n-cycle, where A = 1. Further, with a real R, only a two- fold branching is possible, but a 
complex R allows higher-fold branchings. 

If an n-cycle is not born out of the above branching, we call it an original cycle, e.g., the 4- 
cycle at R' a = 5/4 is born out of bifurcation, see Fig. [jjd), but the other at R a = (3+v / 108)/4 
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is original, see Fig. [TJe), also the discussion after Eq. (7-4). Both types of cycles exist in 
P n (R, A = +1), as separate factors; and the factor responsible for the original cycles, or the 
original factor below, is given by the following formula. 



Theorem 7. The original factor S n (R) of P n (R,X) at the onset is given by 

Pn(RA) 



S n (R) 



n 



cd=n, c>l 



(15) 



n (M =i^(^wc 

where the inner product on the denominator is carried over k from 1 to c that are coprime 
to c. 

We illustrate Theorem [7] through a few examples before giving a proof. For n — 1, 
Si(R) = Pi(R, 1) = —AR — 1 as the denominator is unity. 

For n = 2, P 2 (R, 1) = 4/2-3. But P^R, -1) = -AR+3. So S 2 {R) = P 2 (R, l)/Pi(R, -1) = 
—1. This means that there is no original 2-cycle and the only 2-cycle comes from period 
doubling. 

For n = 3, P3(R, 1) = — (AR — 7)(16i? 2 — AR + 7), whose second factor is equal to 
(-AR + ±=f^) (-AR + ±±^&) = Uk=i,2 p i( R , e 2km / 3 ). Thus S 3 {R) = -AR + 7. 

For n = 4, P 4 ( J R, 1) = (4R- 5)(16i? 2 + 8i? + 5) [(4fl- 3) 3 - 108] . But P 2 (R,-1) = 4i?-5 
(for c = 2) and n*=i,3 p i(^> e ^ i/2 ) = (~ 4/? ~ 2i ~ 1)(— 4i2 + 2z - 1) = 16i? 2 + 8# + 5 (for 
c = 4). Dividing P^(R, 1) by the two factors yields S^(R) = (AR — 3) 3 — 108, whose only 
real root R— (3+ v / 108)/4 corresponds to the onset of the original cycle. Note R = 5/4 is 
excluded from S^R) as it comes from period-doubling the 2-cycle. 

We now prove Theorem [7j Suppose n = cd, we have, from Eq. (J3j) , 

3^2+1 — Xd+l+1 = ~{ x l + x d+l){xi — Xd+l)- 

We apply the equation to I = 1, . . . , m, and the product is 



x m+l ~ x m+d+l — ( — 1) 



1=1 



(X\ — Xd+l)- 



We now set m to 0, d, . . . , (c — 1) d in this equation, add them together, eliminate Xi — a^+i 
(which is nonzero in a cycle), and 

c-1 c'd 

£(-1)^11^ + *<w) = 0. (16) 



c'=0 



Note Eq. (16) holds for every divisor c of n (c > 1). We now have 
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Theorem 8. An n-cycle and a shorter d-cycle (d\n, d < n) intersect only at the onset of 
the n-cycle, and Ylt=i f'( x k) — (— 2) d Xi ...x^ is a primitive (n/d)th root of unity there. 

Proof. At the intersection of the n- and (i-cycles, xi repeats itself after d steps, so x^+i = xf, 



and Eq. (16) becomes, 

l + q + --- + q c ~ l = 0, (17) 
where q = (—2) d xi . . . Xd- Multiplying Eq. (17) by q — 1 yields 1 = q c = (— 2) n x\ . . . x n . So 



the n-cycle is at its onset. 

Further q is a primitive cth root of unity. Suppose the contrary: q = e 2km ^ c and (k, c) = 
g > 1, then by C\ = c/g, k\ = k/g, we have 

g = e 2faffi/ci . (18) 



Similar to Eq. (17), we can apply Eq. (16) with c — )• g and d — > dc%, and 



1 + qi + ■ ■ ■ + qi 9 - 1 = 0, 



where q 1 = (-2) dci x x ... x dci = q Cl . But by Eq. (18), q Cl = e 2kl7ri = 1, and 1 + qi H h 



qi 9 1 = g > 0; a contradiction. □ 

Remark 1. The only real q is q = — 1 for c = 2, i.e., a period-doubling. On the complex 
domain, however, we can have a c-fold branching with c > 2, which corresponds to a contact 
points between "bulbs" in the Mandelbrot set, see Fig. [2j 

By Theorem [8j P n (R, A) at the onset point includes Pd(R, e 2fc7r4//c A 1 / c ) for every possible 
combination of k and c, such that c) = 1, c|n, and c > 1. Dividing the factors from 
P n (R, A) yields Theorem [7| 

I. Degrees in R 

Theorem 9. The degrees in R of A n (R,X), P n (R, A) and S n (R) are 

deg R A n (R, X) = J2 <t>{n/d)2 d -\ (19a) 

d\n 

deg R P n (R, A) = ^u(n/rf)2 d - 1 = /3(n), (19b) 

c(|n 

deg^ n (i?) = /3(n)- ^ (3(d)<f)(n/d). (19c) 

d|n,d<n 
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Proof. We first prove Eq. ( 19 1). We recall the subscript p of C p denotes a sequence of indices 
k in the generating monomial \\ k x^ k ■ But for a square-free cyclic polynomial, each k occurs 
no more than once, so p also represents a set of indices, e.g., p = 1 represents {1}, {C\ = 

x\ + h x n , generator: x\) and p = 13 represents {1, 3} (C13 = £1 £3 + X2 £4 + V x n X2, 

generator: X1X3, assuming n > 4); more examples are listed in Table [T]). In this proof, we 
shall also use p to denote the corresponding index set, \p\ the set size, i.e., the number of 
indices in the set, and p = {1, . . . ,n}\p the complementary set. Obviously, \p\ + \p\ = n. 
Further, we will include p that correspond to alternative generators of the same cyclic 
polynomial, e.g., we allow p = {2}, {3}, . . . , {n}, although they represent the same cyclic 
polynomial C\ as p = {1}. 

Next, we recall the matrix elements T pq (R) arise from the square-free reduction of 
A(x)C p (x) = (-2) C p (x). A single replacement x\? — )• R — Xk+i produces two new 

terms: in the first, x\? — > R, and in the second, x^ 2 — » —Xk+i- We call the two type 1 and 
type 2 replacements, respectively. If a monomial term t(R,x) results from li type 1 and I2 
type 2 replacements during the reduction of a term s(x) in A(x)C p (x), then the degrees in 
x, for any Xk, of s(x) and t(R, x) are related as 

deg x s(x) - deg x t{R, x) = 2l x + l 2 . (20) 

Similarly, the degrees in R satisfy 

de Si? s ( x ) - deg i? t( J R,x) = 

but since deg R s(x) = 0, 

deg Ji t( J R J x) = Z 1 . (21) 

Now if the monomial t(R,x) settles in the qth column of the matrix T(R), as part of 
T pq (R)C g (x) in Eq. ([5]), then t(R,x) must be a generator of C 9 (x); so 

deg x t{R,x) = \q\. (22) 

Since s(x) is part of A(x)C p (x), we have 

de gx s ( x ) = de g x A (x) + deg x C p (x) = n + \p\. (23) 



From Eqs. (|20j), (J22J), (J23J), we get 

n+\p\- \q\ = 2k+l 2 , 
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and 



h = (n+ \p\ - \q\ - l 2 )/2 <(n+ \p\ - \q\)/2. 



(24) 



By Eq. (21), we get 



deg R T pq (R) = max{deg R t( J R,x)} = max^J < (n + \p\ - \q\)/2, 



j pi 



where the equality holds when all replacements are type 1 (Z 2 = 0). 

Finally, each term of the determinant A n (R, A) = | A I— T(i?) | is given by (— l) s JT [A 5 P 
Tpq(R)\, where p runs through rows of the matrix and {q} is a permutation of {p}, with 
(— l) s being the proper sign. Summing over rows under this condition yields 

deg R A n (R, A) = max^degT^GR) < N(n) (n + \p\ - \q\)/2 = nN(n)/2, 



where equality can be achieved if q = p in every row. By Eq. (11) we have Eq. (|19p,). The 
first few values are 1, 3, 6, 12, 20, 42, 70, 144, 270, 540, 1034, 2112, 4108, . . . , starting from 
n = 1. 



To show Eq. (19d), we take the degree in R of Eq. (14). So ^2 d i n P{d) (n/d) = N(n)n/2, 
whose inversion is f3(n)/n = Y2d\n ^{ n /d) N(d)/2 = L(n)/2. The last step follows from 



inverting Eq. (13). The first few values are 1, 1, 3, 6, 15, 27, 63, 120, 252, 495, 1023, 2010, 



4095, . . . , starting from n = 1. 



Eq. (19") follows directly from taking the degree in R of Eq. (15). The first few values 
are 1, 0, 1, 3, 11, 20, 57, 108, 240, 472, 1013, 1959, 4083, . . . , starting from n = I. □ 



Eqs. (19d) was long known [16], and ([T^) was recently derived |17j . 



III. HENON MAP 



We now extend the method to the Henon map |18j : 



x k+1 = 1 + y k - a x k 2 , y k+1 = bx k . 



(25) 



We change variable x k <— ax k , y k <r- ay k , and 



Xk+i = a + y k - x k 2 , y k+1 =bx k . 



(HI) 
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Since neither a nor b is changed during the transformation, Eq. (25) and Eq. (25 ) share the 
same onset and bifurcation points in terms of a and b. We also see that if b — > and a — > R, 
Eq. (25 ) is reduced to the logistic map Eq. rt3j). 

Since y k = bxt-i, we can ignore y k and work with cyclic polynomials of x k only, the 
square free reduction is now Xk 2 — > a + bxk-i — x k+ \. 



The stability of Eq. (25 ) can be found from the Jacobian matrix 



Jb(Xk) 



dxk+i/dxk dxk+i/dy k 
dy k+ i/dx k dy k+ i/dy k 

The eigenvalue A of the composite Jacobian Jb{x\) ■ ■ ■ Jb{ be computed from 



H 


' -2x k l\ 




V b °/ 



A I - J b {xi) ■ ■ ■ J b {x n ) I = A 2 - 0(x) A + {-b) n = 0, 



(26) 



where 9(x) and (— b) n are the trace and determinant of the matrix product Jb(x\) • • • Jb(x n ), 
respectively |19j . In a stable cycle, the magnitude of A cannot exceed 1; so we replace A by 



+1 or —1 in Eq. (26) to obtain the onset or bifurcation point, respectively. Eq. (26) is the 
counterpart of Eq. ^ . 

Since cyclically rotating matrices in a product does not alter the trace, G(x) is a cyclic 
polynomial of x = {x k }. Thus, we can use 0(x) to list Eqs. ^ and then replace O(x) by 
A + (-b) n /X or ±[l + (-&)"] in Eq. Q to complete the solution. 

We computed the polynomials of a and b at the onset and bifurcation points for n up 
to 9. The polynomials of a and b at the onset and bifurcation points, as well as 0(x), are 



listed in Table VI for small n (for larger n, see the website in Section V). 



IV. CUBIC MAP 

We now study the following cubic map [1] 

x k+ i = f{x k ) =rx k - x k 3 . (27) 

Since the new replacement rule 

x k 3 -> r x k - x k+1 (28) 

no longer eliminates squares, we must extend the basis set of cyclic polynomials from the 
square-free ones to the cube- free ones, in using Eq. We include in the basis of expansion 
Cm = xi 2 x 2 + x 2 2 x 3 H h x n 2 xi (n > 3), but not Cm 2 = Xi 3 x 2 + x 2 3 x 3 H + x n 3 Xi. 
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TABLE VI. Onset and bifurcation polynomials of the n-cycles of the Henon map. 



e(x) 



Onset P„(a,6,+1) 



Bifurcation P n (a, 6, — 1) t>i 



-2x, 



4Ci2 + 2b 



-A-{b-l) 2 
A-3(6- l) 2 



-A + 3(6-1) 



A — B 



(2) 
5,-6 



-8C123 — 2bC\ 



-A + B. 



(2) 
7 : 10 



A - ±B (2) 



+ f (b 2 -iy 



A 6 + 2Bf\A 2 - 9B^L 6 - 7 A 



(4) 



•95, 



(6) 
,6,2,-10 



4 16C7i234 + 46Ci2 + 26 2 



~32Ci2345 



-86C123 - 2b 2 d 



( 



A - B, 



(2) 



A 6 - 4A 5 B$ + A*B W 



47,68,38 



{[A+(6+l) 2 ] 2 +4(6 2 -l) 2 } 
{ [A - 3(6 + l) 2 ] 3 - 108(6 - 1) 2 (6 + l) 4 } 



0430(6) 

O^i ^47 62,-83,-212 



r(4) 

- D 17,12,-6 



(4) 



1 R (6) a , r(4) r (4) 

1" -°9, -14, 70, -80^ "r - D l,l,l- D 31,-89,121 



(4) 



A 11 - B ( £ U A 10 + 2B. 



208,206,377 



.4" 



757287226- 28629151 ) 



17,12,-6^31,-60,-186 



./|15 1 Or(2) 4 14 _ d(4) 4 13 

A -f £D XQ X J\ -° 448, 60, 669^ 



1291467969 



512C1...9 

26 4 Ci 



(A 6 + 126A 5 



A 240 - sb^a 239 



)(A 6 



ABf\A 5 



? (4) 4 2 

-•29525, 1019, 58188^ 



120670698649 . . . 712084645033 



2\ a 251 



(A 252 - 504(1 + b 2 )A 



~r ^^31500, 16,62067^ 



250 



-5842146539 . . . 9260477441 



t Definitions: A = 4a, B$ = p(b 2 + 1) + qb, B^% r = p(b 4 + 1) + q(b 3 + 6) + r6 2 , B^% r>a = p(b 6 + 1) + <j(6 5 + 
6)+r(6 4 + 6 2 ) + s6 3 , .... 

* The onset polynomials for n from 1 to 4, and the bifurcation polynomials for n from 1 to 3, agree with 
those in ref. |19j . 



However, we only need the cube-free cyclic polynomials of even degrees in x to solve the 



problem, because Eq. (27) contains only linear and cubic terms, a cyclic polynomial with an 



odd (even) degree in x can never be reduced to one with an even (odd) degree by Eq. (28). 
For technical reasons, we will not use polynomials of odd degrees, because the map allows 
a symmetric 2n-cycle: x%, x 2 , ■ ■ ■ , x n , —xi, — x 2 , ■ ■ ■ , — x n (see Fig. [3J, which makes all odd 
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cyclic polynomials zero, e.g., C\ = x\ + %2 + ■ ■ ■ + x n — X\ — £2 — ■ • • — x n = 0. Thus, the 
zero determinant condition, similar to that in Eq. ([7]), would be useless for these cycles, if 
the odd-degree polynomials were used. 



(a) 



n= 1, r=2.2 



(b) 



n = 2, r= 2.8308 



(c) 



n = 3, r= 2.7003 




FIG. 3. Odd-cycles of the cubic map. 



We therefore have a theorem similar to Theorem [T] 



Theorem 10. For the cubic map Eq. (21), any cyclic polynomial K(x) of an n-cycle orbit 
x = {xi, . . . , x n } with an even degree in x is a linear combination of the even cube-free cyclic 
polynomials C p (x); 

ff(x) = £/„(fl)C p (x), 

where B = {0, 11, 12, 13, ... , 1122, 1123, . . .} is the set of indices of all even cube-free cyclic 
polynomials, and f p (R) are polynomials of R. 

With the above change, the rest derivation is similar to that of the logistic map. The 
new A(x) should be YYk=i f'( x k) = llfc=i( r — 3x fc 2 ). The polynomials of r at the onset 



and bifurcation points for some small n are shown in Table VII (general A) and Table VIII 



(A = ±1); for larger n up to 8, we have saved the data on the website in Section [VJ The 
representative r values are listed in Table IX For complex r and x, we have plotted Fig. [4] 
for regions of stability. 



A. Counting cycles 



We now compute the number of the even cube-free cyclic polynomials by establishing a 
one-to-one correspondence between the cube-free cyclic polynomials and the ternary neck- 

26 



TABLE VII. Minimal polynomials P n (r, A) of the n-cycles of the cubic map Eq. (27). 



n P n (r, A) t 

1 (A-r)(A + 2r-3) 

2 [\^r^Sf\ (A + 2r 2 - 9) 

A 4 + 2(r + 6)(r 2 - 9)A 3 - 6(8r 6 + 12r 5 - 66r 4 - 81r 3 + 54r 2 - 243r - 729)A 2 + 2(r 2 - 9)(l6r 7 - 252r 5 

3 

-216r 4 + 648r 3 + 972r 2 + 2187r + 4374)A + 2r(2r - 3)(2r + 3) 2 (r 2 - 9)(2r 2 - 9) 2 (4r 2 + 9) + 531441 
[(A -8r^J^^ 1 [A 8 + 2(5r 4 - 324)A 7 + 2(112r 8 - 1296r 6 + 3807r 4 - 91854)A 6 

18075490334784r 8 - 61004779879896r 6 - 411782264189298r 4 + 18530201888518411 



t The factors from odd-cycles (see Section 



IV B 



are struck out. 



$ The n = 4 odd-cycles satisfy (A odd + 9) 2 - 2(A odd - 27)r 2 - 8r 4 = 0, where A odd = ±V\. 



laces, in which each bead of the string is assigned a number 0, 1, or 2, instead of just 
or 1. For example, the necklace 212001- •• corresponds to Cn2336(x), whose generator is 
X\ 2 X2X3 2 Xq: the first bead is 2 for x\ 2 , the second is 1 for X2 1 , the third is 2 for X3 2 , and 
the sixth is 1 for xq 1 . A necklace is even, if the corresponding cyclic polynomial has an even 
degree in x. This means that the sum of numbers (0, 1, or 2) on the beads of the necklace, 
which equals the degree in x of the polynomial, is also even. 

Theorem 11. The number of the even ternary necklaces or the cube-free cyclic polynomials 
for the cubic map of even degrees in x is given by 

3 d - 1" 



N e (n) = -J2 <M 



n 

cd=n 



odd(c) 



(29) 



where odd(c) = [l — (— l) c ]/2 is 1 if c is odd or if even. 

Proof. We first show that the number of even ternary strings is (3 n + l)/2. Consider the 
generating function 

n 

Z({x 1 ,X 2 , ■ ■ -,X n }) = + Xi + X 2 ), 



i=l 

.2 



where 1, Xi, and x~ correspond to the states that bead i taking the number 0, 1, and 2, 
respectively; and the product over n sums over states of independent beads. In the expansion 
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TABLE VIII. Onset and bifurcation polynomials of the n-cycles of the cubic map Eq. (27). 



Onset P n (r,+1) 1 



Bifurcation P n (r, — 1) t 



-(r-1) 2 
-2(r-2)(r + l)(r + 2) 2 

(r 2 + r + l)(4r 2 - 14r + 13) 

(4r 8 + 16r 7 - 35r 6 - 206r 5 - 113r 4 

+ 376r 3 + 715r 2 + 1690r + 2197) 

-8(r 2 - 8)(r 2 - 5)(r 2 + l)(2r 2 - 6r + 5) 
(2r 2 + 6r + 5)(2r 4 - 13r 2 - 25) 2 
(1024r 22 - 32512r 20 + 402304r 18 



-2364832r 16 + 5389924r 14 + 9715769r 12 



-73067038r lu + 58934785r 8 + 235761152r b 
-160907264r 4 - 671088640r 2 - 2097152000) 



-(r + l)(r-2) 
-(r 2 -5)(2r 2 +6r + 5) 



16r 12 + 24r n - 288r 10 - 434r 9 



+1539r 8 + 2358r 7 - 1434r 6 - 2556r 5 - 8541r 4 
-11816r 3 + 15288r 2 + 24696r + 38416 

- (16r 8 - 216r 6 + 410r 4 + 2142r 2 + 1681) 
(8192r 32 - 387072r 30 + 7834624r 28 - 88 03 1 232r 26 
+ 585876512r 24 - 2158227720r 22 + 221136 1 3 1 2r 20 



15958823175r 18 - 68871388441r 16 + 59290039854r 14 



234882618673r 12 - 524807876277r 10 - 72612143404r 8 



308406843576r 6 + 1539579145957r 4 - 7984925229121) 



- 8192(r 4 + l)(2r 4 - 18r 2 + 41) 
(8r 4 - 48r 3 + 108r 2 - 108r + 41) 
(8r 4 + 48r 3 + 108r 2 + 108r + 41) . . . 
(22016722240 . . . 14954924752896r 3108 



-(17592186044416r 80 + • • • + 
144564714832407908402064153121600801) 
(4841528421712030 . . . 03048551060078592 
25263420710 . . . 173884723232001) 



,3200 



180097954 . . . 2522413056 x 10 4UU ) 

t Although P n (r, A) is the minimal polynomial for a general A, it may contain a pre-factor (e.g., —2 in the 
n = 2, A = +1 case). 



of Z(x), each term (which takes the form xi ei X2 S2 ■ ■ -x n en , with = 0,1,2) represents a 
unique ternary string eie2 . . . e n , which is even, if e\ + ti + • • • + e n is so. By setting 
%i = X2 = ■ ■ ■ = x n = 1, Z = 3 n equals the total number of ternary strings. By setting 
x\ = X2 = ■ ■ ■ = x n = —1, a term representing an even (odd) ternary string is +1 (—1); and 
Z — (1 — 1 + l) n = 1 equals the difference between the number of even strings and that of 
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TABLE IX. Smallest positive r at the onset and bifurcation points of the n-cycles of the cubic 
map. 





Onset* 


Bifurcation* 


#* 


rat 


Onset* 


Bifurcation* 


#* 


1 


l.OOOOOOOOOOi 


2.0000000000i 


2 


6 


2.3334877526 30 4 


2.3355337580 336 


56 


2' 


2.0000000000i 


2.2360679775 2 


2 


6' 


2.4608286739i2 


2.4657090579 336 


4 


3 


2.4504409645 8 


2.4608286739i2 


4 


7 


2.3729872678io8o 


2.3732727868io92 


156 


4 


2.5478350393 22 


2.5488312193 32 


8 


8 


2.3525990555 3 io8 


2.3527637793 320 o 


400 


4" 


2. 2360679775 2 


2.2880317545 32 


2 


8' 


2.5488312193 32 


2.5493247379 3 2oo 


8 


5 


2.3939250274n2 


2.3957922744i2o 


24 


8"' 


2.2880317545 32 


2.29 9 2 2 793 9 7 3 2oo 


2 



' ', ", or '" means a cycle under the first, second, or third successive period-doubling, respectively. 

* The subscripts are the degrees of the corresponding minimal polynomial. 

* The number of similar cycles (for n > 1, only half of them have positive r). 



odd strings. Thus, the average (3 n + l)/2 gives the number of even ternary strings. 



TABLE X. Two ways of counting the N e (2) = 4 ternary necklaces for n = 2. 





m = 1 


m = 2 


c = n/d 


T(d,c) * 


0(c) 


<j>{n/d)T{d, c) 


d = 1 


0x2 1x2 2 X 2 




2 


3 


1 


3 


d = 2 


00 11 22 


20(02) * 


1 


3 + 1-2 = 5 


1 


5 


TO Em| £ M|n<M§) 


1- (1 + 1) = 2 


2-1 = 2 




1 


N e (n)-n 


3-2 


1 • 2 




4-2 = 3 + 5 



Bold strings are necklaces; others are their cyclic versions. 

The subscript of a string means the number of repeats; e.g., 1 X 2 means 1 repeated twice, or 11; 

* T(d, c) ee 3 d - odd(c)(3 rf - l)/2, which is 3 d if c is even, or (3 d + l)/2 if c is odd. 

* The odd binary strings 10, 01, 21, and 12 do not contribute to the sum, and are excluded from T(d, c). 



The rest counting process is similar to that in Section II G we construct the sum nN e (n) 
by two ways, as exemplified in Table |Xj In the first way, for a fixed d (d\n), if c = n/d 
is even, we count all ternary strings whose period m divide d, but if c is odd, we count 
only ternary strings whose first d beads are even (because repeating an odd string an odd 
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1 - 



- 




-2 



-1 





Re(r) 



FIG. 4. Stable regions of the n-cycles of the cubic map Eq. (27) with a complex R; obtained from 
p n (U, e i(f> ) = with <j> € [0, 2vr], cf. Fig. § 



c . 



number times does not yield an even string); in either case, we multiple the result by 
The process for a fixed d corresponds to a row of Table |X} Repeating the process for all d 
gives J2cd=n T(d, c) 0(c), where T(d,c) is the total number 3 d of period-fi ternary strings if 
c is even, or the number (3 d + l)/2 of even strings if c is odd. 

In the second way, we look at the contribution from each necklace to the above sum. An 
even period-m necklace contributes a total of m x J2 m \ dd \ n 4>{n/d) = n [the multiplier m is 



for the m cyclic versions, cf. Eq. (10)], while an odd necklace contributes nothing. Thus, 
the sum equals N e (n) ■ n. The process for a fixed necklace corresponds to a column of Table 

m 



So 



N e (n) ■ n = W 



cd=n 



which is Eq. (30) after we divide both sides by n. 



□ 



N e (n) = 2, 4, 6, 14, 26, 68, 158, 424, . . . , starting from n = 1. 

The characteristic polynomial A n (r, A) from the determinant equation has a degree N e {n) 
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in A. Again, it encompasses the factors for the n-cycles and the shorter o?-cycles, as long as 
d\n. The minimal polynomial for the n-cycles can be obtained by Theorem [3] with proper 
substitutions [R — > r, A n (R, A) — > A n (r, A), etc.] . The degree of the polynomial is given by 



Theorem 12. The degree in A of the minimal polynomial P n (r, A) of the n-cycles is 

3 d -r 



LJn) 



n 

cd=n 



1 + odd(c) 



2 

where odd(n) = [1 — (— l) n ]/2 is 1 for an odd n, but for an even n. 



(30) 



Eq. (30) follows from the inversion L e (n) = Yld\n l J '( n /d) N e (n), after some algebra, as 
shown in Appendix [Bj L e (n) = 2, 2, 4, 10, 24, 60, 156, 410, . . . , starting from n — 1. This 
is also the number of the n-cycles [13] . Note that, for n > 1, half of the cycles have negative 
r, and the Xk are imaginary. However, in a transformed map, 

= r z k (1 — ^fc 2 ), 



which differs from Eq. (27) by a;*; = y/rzk, Zk in the negative-r cycles are real. 

Following a similar proof to Theorem [9j we find the corresponding degrees in r of the char- 
acteristic polynomial A n (r, A) and minimal polynomial P n (r, A) of the n-cycles are nN e (n) 
and nL e (n), respectively. 

B. Odd-cycles 

Because of the symmetry f(—x) = —f(x), the minimal polynomial P n (r, A) is subject 
to factorization for an even n. If X(„/2)+i = — X\, then xi, . . . ,x n /2, —xi, . . . , —x n /2 is an 
n-cycle, for x n+ i = — X( n /2)+i = X\. We call such a cycle an odd-cycle, see Fig. [3] for 
examples. Odd-cycles satisfy a polynomial of lower degrees in A, which causes the factor- 
ization. Suppose A odd (x) = Y\k=i f'( x k) i n the odd-cycle satisfies P° d ^(r, A odd ) = [where 
A odd is the value of A odd (x), and A odd = ±VX], then P°f£(r, y/\)Pffi(r, — \/X) is a factor 
of P n (r, A). For example, by solving the n = 2 odd-cycle, we have — Xi = rxi — X\ 3 , or 
xi 2 = r + 1. Since A odd = r - 3xi 2 , P 1 odd (r, A odd ) = A odd + 2r + 3. Now the factor for 2-cycles 
(see Table |Vn|) is P 2 (r, A) = - [(2r + 3) 2 - A] (A + 2r 2 - 9), whose first factor is indeed 
P 1 odd (r, >/A) P 1 odd (r, -VA) = (VX + 2r + 3)(-v / A + 2r + 3). The example for n = 4 is shown 
in Table EU 
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V. SUMMARY AND DISCUSSIONS 



We now summarize the algorithm for a one-dimensional polynomial map. First, we list 
Eqs. (|6j) with A(x) = YYk=i f'( x k)- This step populates elements of the matrix T(r), where r 
is the parameter of the map. The determinant A n {r, A) = I A I — T(r) | , with A being +1 and 
— 1, then gives the characteristic polynomial at onset and bifurcation points, respectively. 
To filter out factors for the shorter d-cycles with d\n, we repeat the process for other divisors 
d of n and then apply 

When implemented on a computer, it is often helpful to evaluate A n (r, A) by Lagrange 
interpolation, that is, we evaluate A n (r, A) at a few different r, e.g., r = 0, ±1, ±2, . . ., then 
piece them together to a polynomial. The strategy also allows a trivial parallelization. 

The algorithm (implemented as a Mathematica program) was quite efficient. For the 
logistic map, the bifurcation point for n = 8 took three seconds to compute on a desktop 
computer (single core, Intel® Dual-Core CPU 2.50GHz). In comparison, the same problem 
took roughly 5.5 hours [I] using Grobner basis and 44 minutes in a later study [12J. To be 
fair, using the latest Magma, computing the Grobner basis took 81 and 14 minutes, on the 
same machine for Eq. Q and Eq. ([3]), respectively; even so, our approach still had a 200-fold 
speed-up. 

The exact polynomials of these maps are generally too large to print on paper, e.g., 
the polynomial for the logistic map with n = 13 takes roughly seven megabytes to write 
down. We therefore save the polynomials and programs of the three maps on the web site 
http: / /logperiod. appspot.com. 



ACKNOWLEDGEMENTS 



I thank Drs. T. Gilbert and Y. Mei for helpful communications. Computing time on 
the Shared University Grid at Rice, funded by NSF under Grant EIA-0216467, is gratefully 
acknowledged. 

32 



Appendix A: Simple derivation of 4-cycles 



The polynomials for the 4-cycles permits a short derivation. We first list the explicit 
equations: 



(Ala) 
(Alb) 
(Ale) 
(Aid) 



X2 


= R- 


x ii 


X3 


= R- 




£4 


= R- 






= R- 


x 2 



[Eq. flAlfr) - Eq. ([Alp)] x [Eq. flAljp) - Eq. (|ATji)l yields 1 + (aci + £3) (£2 + x 4 ) = 0, since 
Xi 7^ x 3 , x 2 7^ x 4 . Hence, with 2/1 = x\ + x 3 , 2/2 = x 2 + x 4 , z = yi + y 2 , we have 



2/12/2 

.2 1 „.2 



2/i 2 + 2/2 = (2/1 + 2/ 2 ) 2 - 2j/ij/ 2 = z 2 + 2 
2/? + 2/2 = (2/1 + 2/2) 3 - 32/i2/ 2 (2/i + y 2 ) = z 3 + 3z. 



(A2a) 
(A2b) 
(A2c) 



Multiplying Eq. (All) by x\ or 23, then summing over cyclic versions yields 



2/12/2 = Rz- [(xl + xl) + (x\ + x\)] , 

2/12/2 = Rz - [xix 3 (xi + x 3 ) + £2^4(^2 + x 4 )] ■ 



(A3a) 
(A3b) 



From Eq. flA3}i) + 3 x Eq. (|A3Jd), we have 4 2/12/2 = 4i?z - (2/? + and by Eqs. flA2| ), 

^ 3 -(4i?- 3)^-4 = 0. 
Since 2xiX 3 = yf — (xf + xl) = y\ — 2i? + 2/2, and 2x 2 x 4 = y\ — 2R + 2/1, 

X = X1X2X3X4 = - z) + (R 2 - R + 1), 



(A4) 



(A5) 



where we have used Eqs. (A2) and (A4) to simplify the result. Dividing the polynomial in 
Eq. (|A§ by that in Eq. (A5} yields z = (R 2 — 3R — X + 1)/(R 2 - R + X - 1), and plugging 
it back to Eq. (gHJ gives i? 6 - 3i? 5 + (3 + X)(i? 4 - R 3 ) + (1 - X)(2 + X)R 2 + (1 - X) 3 = 0, 
which is the same as the first factor of Eq. (gj) with X = A/16. 
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Appendix B: Proof of Theorem 



12 



Here we prove Theorem 12 [or Eq. (30)] for the cubic map. Similar to the logistic map 



case Eq. (13), we have, for the cubic map, 



N e {n) = Y,L e {d). 

d\n 

Thus, we only need to inverse this equation to obtain L e (n). But owing to the complexity 



of Eq. (29), we need the Dirichlet generating function to simplify the result. 



For a series a(n), the Dirichlet generating function is defined as 

oo 



a\n)n 



n=l 



In Table |XI[ we list the generating functions of some common series, and define a few new 
ones for N e (n) and L e (n), etc.. 

TABLE XI. Dirichlet generating functions for ternary necklaces. 



a(n) 




a(n) 


G a (s) = J2 n a(n)/n s 


1 




C(s) f 




N e (n) 


G N (s) i 


S n ,i 




1 t 




nN e (n) 


G N (s-l) t 


H(n) 




COT 1 f 




L e (n) 


Gl(s) { 


4>{n) 




C(*-l)/C(*) 


t 


nL e (n) 


G L (s-l) * 


odd(n) 




C(s)(l-2" s ) 


3" 


t(s) i 


/z(n) odd(n) 




'C(s)(l-2- s ) 


-l 


4>(n) odd(n) 


f(s-l) l-2~ s + 1 
C(s) 1-2- = 



' C( s ) = n s is the zeta function. See ref. [T3] for proofs. 
* The sum is truncated at a large M to avoid divergence. 



The generating function has an important property: G 7 (s) = G ct (s)G j a(s), if and only if 
7(71) = J2d\n a ( n /d) (3(d) [13]. Thus, the terms of the sum ^2 d \ n a{n/d) (3(d) of two sequences 
a and (3 can be readily found from expanding the generating function. 

Another fact is if G a (s) is the generating function of a(n), then G a (s — 1) is that of 
na(n), for a{n)/n s ~ l = [na(ri)] /n s . Thus, the generating function of n is C( s— 1)> an d that 
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of nN e (n) is Gn(s — 1) [((s) is the generating function of 1, and Gn(s) is that of N e (n), see 
Table 



XI 



We now compute the generating function of fi(n)odd(n). First, recall the generating 
function G M (s) of /i(n) is 

fi{n) 



n H 



n=l p 

where p goes through every prime. The follows directly from expanding the product and the 
definition of fi(n), which is —1 to the power of the number of distinct prime factors. The 
same reasoning applies to G ^^(s) = Y^ n odd ^(n)/n s with the only difference being that all 
multiples n of 2 are absent. So 



G>,odd(s) - } j 



n(n) odd(n) 



n=l 



rr 



E 

n odd 



H{n) 



rr 



p>3 



G M (s) ( 1 



1 



-i 



COO i 



-1 -1 



Comparing the two formulas yields 

Gfj, )0 dd{ s 

Similarly, we can compute the generating function of odd(n) as 

n=l n odd x ' 

This can also be derived by taking the generating function of both sides of the identity: 
Tld\n n(d)odd(d) odd(n/d) = 5 n> i [which is a modification of E^ d | n /i(<i) = 5 nj i] . It follows 
that the generating function of nodd(n) is G dd(s — 1). 

The generating function G^m^s) of 0(n)odd(n) can be computed by taking the gener- 
ating function of both sides of the identity 

nodd(n) = (f>(d) odd(d) odd(n/d), 

d\n 

i.e., if n is even, then both sides are 0; if odd, then n = E^ d | n 0(<f). So 

r , , _ G odd (s - 1) _ C(s - 1) 1 - 
* Md[S) ~ G odd (s) ~ COO 1 - 2- ' 
We can now compute the generating function Gn(s) of N e (s). By multiplying n to both 



sides of Eq. (29), and taking the generating function, we find that 

i( 8 )-C(a) C( S -l)[^) + C(^)(l-2- +1 ) 



G N (s - 1) = G^(s)t(s) - G^odd(s)- 



2C(s)(l-2" 
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where the left side Gn{s — 1) is the generating function of N e (n)n, and formulas in Table 
IXJhave been used. 

Finally, we take the generating function of both sides of N e (n) = ^2 d \ n L e {d): 

G N (s) = ((s)G L (s), 



and 



G L {s-l) 



G N (s-l) 



t(s)-C(s) , too- coo ^ / \ 



C(s-l) 2 COO (1- 2-) ^ 2 

Comparing the coefficients of the nth term (n <C M), we find 

3 d - 1 



nL e (n) = 8 n>1 + 22 M c ) odd 



cd=n 



which is Eq. (30) [also note 5 n ^ = J2 c i n ^(c 
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